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ABSTRACT 

We have constrained the age of the globular cluster S312 in the Andromeda galaxy (M31) by comparing 
its multicolor photometry with theoretical stellar population synthesis models. This is both a check 
on the age of this globular cluster, as well a check on our methodology. Main-sequence photometry 
has been the most direct method for determining the age of a star cluster. S312 was observed as part 
of the Beijing- Arizona- Taiwan-Connecticut (BATC) Multicolor Sky Survey from 1995 February to 
2003 December. The photometry of BATC images for S312 was taken with 9 intermediate-band filters 
covering 5000 — 10000A. Combined with photometry in the near-ultraviolet (NUV) of GALEX, broad- 
band UBVR and infrared JHK S of 2MASS, we obtained the accurate spectral energy distributions 
(SEDs) of S312 from 2267 — 20000A. A quantitative comparison to simple stellar population models 
yields an age of 9.5l 'g g Gyr, which is in very good agreement with the previous determination by 
main-sequence photometry. S312 has a mass of 9.8 ± 1.85 x 10 5 M Q , and is a medium-mass globular 
cluster in M31. By analysis of errors of ages determined based on the SED fitting method of this paper, 
secure age constraints are derived with errors of < 3 Gyr for ages younger than 9 Gyr. In fact, the 
theoretical SEDs are not sensitive to the variation of age for ages greater than ~ 10 Gyr. Therefore, 
for globular clusters as old as the majority of the Galactic GCs, our method do not distinguish them 
accurately. We emphasize that our results show that even with multiband photometry spanning NUV 
to K s , our age constraints from SED fitting are distressingly uncertain, which has implications for 
age derivations in extragalactic globular cluster systems. 

Subject headings: galaxies: individual (M31) - galaxies: star clusters - galaxies: stellar content 



1. INTRODUCTION 

Galactic globular clusters (GCs), which are thought 
to be among the oldest stellar objects in the Universe, 
provide vitally important information regarding the min- 
imum age of the Universe and the early formation his- 
tory of our Galaxy. The most direct method for deter- 
mining the age of a star cluster is main-sequence pho- 
tom etry, since the tur n-off is mostly affected by age 
(see iPuzia et al.l [200 2b. and references therein). How- 
ever, this method has only been applied to the Galac- 
tic GCs and glo bular clusters in the sat ellites of the 
Milky Way (e.g., iRich et all f2001h before IBrown et all 
(2004). Generally, extragalactic globular cluster ages 
are inferred from comp osite col o rs an d/or spectroscopy. 
S312, first detected by ISharov I (|1973h ( No.19), then by 
ISargent et ail (|1977f) (No.312=S312) and lBattistini et all 
(|1987ft (No"!379=Bo379). is among the first extragalac- 
tic globular clusters whose age was a ccurately es- 
timat e d by main-sequen ce photometry (jBrown et al.l 
12004) . IBrown et all (|2004h obtained the color -magnitude 
diagram (CMD) below the main-sequence turnoff for 
S312 using the extremely deep images of M31 with the 
Advanced Camera for S urveys (ACS ) on t he Hubble 
Space Telescope (HST). IBrown et all (|2004h estimated 
an age of 10 Gyr fo r S312 by the quantitativ e compari- 
son to isochrones of IVandenBerg etldl (|2006l ). 

This cluster was also observed as part of the galaxy 
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calibration program of the Beijing- Arizona- Taiwan- 
Connecticut (BATC) Multicolor Sky Survey (e.g., 
iFan et alJ 119961 : I Zheng et all Il999f ) in 9 intermediate- 
band filters. Combined with photometry in NUV of 
GALEX, broad-band UBVR and infrared JHK S of the 
Two Micron All Sky Survey (2MASS), we obtained the 
accurate SEDs of S3 12 in 17 filter bands from 2267 - 
20000A. 

Sinc e the p ioneering work of lTinslevI (|1968l . I1972D and 
ISearld (|1973l ). evolutionary population synthesis mod- 
elling has become a powerful tool to interpret inte- 
grated spectrophotometric observations of galaxies and 
their subcomponents, such as star clusters ([Anders et al.l 
l2004f h The evolution of star clusters is usually modelled 
by means of a simple stellar population (SSP) approx- 
imation. An SSP is defined as a single generation of 
coeval stars formed from the same progenitor molecu- 
lar cloud (thus implying a single metallicity) , and gov- 
erned by a given initial mass function (IMF). Glob- 
ular clusters, which are bright and easily identifiable, 
are typically characterized by homo geneous abundance 
and a ge distributions. For example. iBarmbv fc Huchral 
(|2001f ) compared the predicted SSP colors of three stel- 
lar population synthesis models to the intrinsic broad- 
band UBVIRJHK colors of Galactic and M31 GCs, and 
found that the best-fitting models match the clusters' 
SEDs very well indeed. 

In this paper, we constrain the age of S312 by compar- 
ing observational SEDs (Sect. 2) with population syn- 
thesis models in Sect. 3. Our independently-constrained 
result is in ve ry good agreement with the previous de- 
termination of IBrown et all ([200l . We give discussions 
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and summarize our results in Sect. 4. 

2. OPTICAL AND INFRARED OBSERVATIONS OF GC S312 

2.1. Historical observations of GC S312 

S312 was first detected by iSharov I (| 19731) . who 
searched 25 M31 globular clusters and cluster candidates 
using plates taken with the 80/120/240-cm Schmidt tele- 
scope of the Radio Astrophysical Observat ory of the Lat- 
vian A cade my of Sciences. Since th at time. ISar gent et 
(|1977f) and IBattistini et all (|1987f ) confirmed S312 to be 
an M31 globular cluster. 

2.2. Archival images of the BATC sky survey 

Observations of S3 12 were obtained by the BATC 
60/90cm Schmidt telescope located at the XingLong sta- 
tion of the National Astronomical Observatory of China 
(NAOC). This telescope has 15 intermediate-band filters 
covering the optical wavelength range of 3000 — 10000 
A, and is specifically designed to avoid contaminations 
from the brightest and most variable night-sky emission 
lines. Descri ptions of the BAT C photometric system can 
be found in iFan et all (|1996t ). The finding chart of GC 
S312 in the g band (centered on 5795 A) of the BATC 
system with the NAOC 60/90cm Schmidt telescope is 
illustrated in Figure 1. 

We extracted 125 images of M31, taken in 9 BATC 
filters from the BATC survey archive during 1995 
February-2003 December. Table 1 contains the log of 
observations. Multiple images of the same filter were 
combined to improve the image quality. 

2.3. Intermediate-band photometry of S312 

Since S312 is located in the M31 halo, it is easy to do 
the photometry. The intermediate-band magnitudes of 
S312 are determined on the combined images with the 
standard apert ure photometr y, i.e. the PHOT routine 
in DAOPHOT (|StetsorJll987f ). The BATC photometric 
system calibrates the magnitude zero level similar to the 
spectrophotometric AB magnitude system. For the flux 
calibration, the Oke-Gunn primary flux standard stars 
HD 19445, HP 84937, BP +26° 2606, and BD +17°4708 
taken from lOke fc Gunnl (1 19831). were observed during 
photometric nights ( Yan et al.ll2000l ). The results of well- 
calibrated photometry of S312 in 9 filters are summarized 
in the fifth column of Table 1 . 

2.4. GALEX, broad-band and 2MASS photometry of 

S312 

In order to constrain the age of S312 accurately, we 
use as many photometric data points covering as large a 
wavelength range as possible. M31 field was observed as 
part of the Nearly Galaxy Survey (NGS) carried out by 
Galaxy Evolution Explorer (GALEX) in two UV bands: 
far-ultraviole t (FUV) a nd near-ultra violet (NUV) (see 
details from iRev et all 120051 . l2006f) . iRev et alj (j2006l ) 
presented the photometric data for 485 and 273 M31 
GCs in GALEX NUV and FU V, respectively. S312 
was not detected in FUV from I Rev et aD (|2006l . also 
S.-C. Rey, priv. comm.). Using the CCD imaging of 
0.9m telescope at the Kitt P eak National Observatory, 
iReed. Harris k. Harris! (|1994f ) presented CCD BVR in- 
tegrated magnitudes and color indices for 41 globular 
clusters and cluster candidates in the outer halo of M31 



including S312. iBattistini et~aTl (|1987l ) presented UBVr 
photometry from the photographic plates for their most 
sample globular clusters and candidates including S312, 
however they did not give photometric uncertainties. In 
this paper, for S312, we adopted CCD BVR photome- 
try of lReed. Harris fc Harris! (119941). a nd photographic U 
photometry of IBattistini et al.l (|1987h with a photomet- 
ric un certainty of 0.08 mag as suggested bv lGalleti et al.l 

w. 

As pointed out by IWorthevi (|1994D . the age- 
metallicity degeneracy in optical broad-band colors is 
Aage/AZ ~ 3/2, implying that the composite spectrum 
of an old stellar population is indistinguishable from 
that of a younger bu t more metal-rich popu lat ion (and 
vice v ersa) (also see iMacArthur et all l2004f) . Ide Jond 
(|1996| ) showed that this degeneracy can be partially bro- 
ken by addi ng 2MASS infrared p hotometry to the op- 
tical colors. ICardiel et alj (|2003f ) found that inclusion 
of an infrared band can improve the predictive power 
of the stellar population diagnosti cs by ~ 30 times over 
using optical photometry alone. IWu et all (2005) also 
showed that the use of near-infrared colors can break 
the age-metallicity degen erac y. At the s a me tim e, as 
iKissler-Patig et all (|2002f ) and lPuzia etaH (|2002af) indi- 
cated that near-infrared photometry is less sensitive to 
interstellar extinction with respect to the classical opti- 
cal bands, it provides useful complementary information 
that can hel p to disentangle the age-metallicity degener- 
acy (also see lGalleti et all 12004] ) . 

In this paper, we add JHK S photometry from the 
2MASS Point Source Catalogue (PSC) to broad- and 
intermediate-band optical, and ultraviolet photometry 
to constrain the age of S312 to the highest possible 
accuracy. 2MASS presents J, H and K s complete ho- 
mogeneous photometry of 99.998 % of the sky down 
to K a ~ 15.5 (completeness > 99 % for K a < 14.3). 
The near-ultraviolet, broad-band and 2MASS photomet- 
ric data points of S312 are listed in Table 2. 

3. STELLAR POPULATION OF S312 

3.1. Stellar populations and synthetic photometry 

As a check on the age of S312, and also on a check 
of our methodology, we compare its spectral energy dis- 
tributions with theoretical stellar p opulation synthesis 
model s. We used the SSP models of iBruzual fc Charlotl 
(2003) (hereafter BC03) , which have been upgraded from 
the IBruzual fc Charlotl (|1993l . [l996) version, and now 
provide the evolution of the spectra and photometric 
properties for a wider range of stellar metallicities. BC03 
provide 26 SSP models (both of high and low resolution) 
using the 1994 Padova ev olutionary track s, 13 of which 
were computed using the IChabrierl (pOOl ) IMF assum- 
ing lower and upper mass cut-offs of tol = 0.1 Mq and 
mu = 100 M(T) , respect i vely. The other 13 were com- 
puted using the lSalpeteJ (|1955l ) IMF with the same mass 
cut-offs. In addition, BC03 provide 26 SSP models us- 
ing the 2000 Pado v a evo lutionary tracks. However, as 
IBruzual fc Chariot) (|2003l ) pointed out, the 2000 Padova 
models, which include more recent input physics than 
the 1994 models, tend to produce worse agreement with 
observed galaxy colors. These SSP models contain 221 
spectra describing the spectral evolution of SSPs from 
1.0 x 10 5 yr to 20 Gyr. The evolving spectra include 
the contribution of the stellar component at wavelengths 
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Fig. 1.— An imago of S312 in the BATC g band of the NAOC 60/90cm Schmidt telescope. S312 is circled. The field of view of the 
image is 11' X 11'. 



from 9lA to 160/im. In this paper, we adopt the high- 
resolution SSP models co mputed using t he 1994 Padova 
evolutionary tracks and a lSalpeterl (|1955h IMF 4 . We note 
that although the current best constraint of the age of 
the Universe is of order 13.7 Gyr, the SSP models and 
the stellar evolutionary tracks that form their basis have 
been calculated for ages up to 20 Gyr. It is not straight- 
forward to correct for this discrepancy; one would need 
to recalculate all stellar evolutionary tracks for all metal- 
licities. 

Since our observational data are integrated luminosites 
through our set of filters, we convolved the BC03 SSP 
SEDs with the BATC intermediate-, GALEX NUV-, 
broad-band UBVR and 2MASS filter response curves 
to obtain synthetic ultraviolet, optical and near-infrared 
photometry for comparison. The synthetic ith filter mag- 
nitude can be computed as 



-2.5 log 



/ A F A ^(A)dA 
f x ipi(\)d\ 



48.60, 



(1) 



where Fa is the theoretical SED and (pi the response 
curve of the ith filter of the BATC, GALEX, UBVR 

4 We note that because of the slow SED evolution of SSPs at 
ages in excess of a few Gyr, all of the most commonly used spectral 
synthesis models agree very well at these ages. Therefore, the 
choice of IMF is only important for estimating the photometric 
mass of the cluster, and does not affect the determination of the 
age of S312. 



and 2MASS photometric systems. Here, Fx varies with 
age and metallicity. 

3.2. Reddening and metallicity of S312 

To obtain the intrinsic SEDs of S312, t he ph otome- 
try must be dereddened. iBarmbv et al.1 (|2000D deter- 
mined the reddening for each individual cluster using 
correlations between optical and infrared colours and 
metallicity and by defining various "reddening-free" pa- 
rameters using their lar g e dat abase of multi-colour pho- 
tometry. IBarmbv et~aTl (|2000D found that the M31 and 
Galactic GC extinction laws, and the M31 and Galactic 
GC colour-metallicity relations are similar to each other. 
They then estimated the reddening to M31 objects with 
spectroscopic data using the relation between intrinsic 
optical colours and metallicity for Galactic clusters. For 
objects without spectroscopic data, they used the re- 
lationships between the reddening-free parameters and 
certain intrins i c colo urs based on the Galactic GC data. 
IBarmbv et al.l (|2000f ) compared their results with those 
in the literature and confirmed that their estimated red- 
dening values are reasonable, and quantitatively consis- 
tent with previous determi nations for GCs acros s the en- 
tire M31 disc. In particular. lBarmbv et al.l (|2000f ) showed 
that the distribution of reddening values as a function 
of position appears reasonable in that the objects with 
the smallest reddening are spread across the disk and 
halo, while the objects with the largest reddening are 
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concentrated in the galacti c disk. Therefo r e, we adopted 
the reddening values from iBarmbv et al.1 (|2000L also P. 
B armby, priv . comm .) for S312 to be E(B - V) = 0.10 
as iRev et all (|2006T > did. The values of extinction coef- 
ficient R\ are obt ained by interpolating the interstellar 
extinction curve of ICardelli et al.l (| 1989( 1. 

The SEDs for clusters are significantly affected by 
the adopted metallicity, especially for old clusters. So, 
only if the metallicity is kno wn, the age o f a clu ster 
can be constrained accurately. iHuchra et al.l (|1991f ) de- 
rived metallicities for 150 M31 globular clusters includ- 
ing S312, using the strengths of six absorption features 
in the cluster in tegrated spec t ra. The metallicity of 
S312, obtained bvlHuchra et al.l (fl99lh. is [Fe/H] = -0.7. 
iHolland et all (jl997D used HST WFPC2 photometry to 
construct deep color-magnitude diagrams for S312, and 
the shape of the red giant branch gave an iron abun- 
dance of [Fe/H] = —0.53. In this paper, we ado pt for 
comparison [Fe/H] = -0.53 as lBrown et all ([2001 did. 

3.3. Comparison of the age scale 

Before constraining the age of S312, we compare the 
age scale between t he Padova 1994 evolutionary tracks 
(|Bertelli et al.lll994l ) used by BC03 and the 2006 van- 
denBerg isochrone s (|VandenBerg et ail 120061 ) used in 

1 Brown et alj (|2004f ). As an example, we only draw the 
isochrones with 10 Gyr and the solar metallicity. Figure 

2 shows the isochrones between effective tcmperaure and 
bolometric magnitudes. It is clear that the matching is 
very good in the main sequence (MS) and the subgiant 
branch (SGB). As we know, the main-sequence photom- 
etry focuses on the upper end of the MS and SGB, since 
uncertainties in models of stars on the red g iant branch 
(RGB) are large. So, if lBrown et all (|2004( l had fitted 
the color-magnitude diagram (CDM) of S312 with the 
Padova 1994 evolutionary tracks, its age would have been 
close to 10 Gyr. 

3.4. Fit results 

We use a x 2 minimization test to examine which 
BC03 SSP models are most compatible with the observed 
SEDs, following 

X =2^ -s , (2) 

i=i 1 

where mJf od {t) is the integrated magnitude in the ?th 
filter of a theoretical SSP at age t, m 1 " tr presents the 
intrinsic integrated magnitude in the same filter and 

°i = ^obs,! + °mod,i- (3) 

Here, a^ hs i is the observational uncertainty, and cr^ lod i is 
the uncer tainty associated wit h the model itself, for the 
ith filter. I Chariot et alj (|1996l ) estimated the uncertainty 
associated with the term cr^ od i by comparing the colors 
obtained from different stel lar evolutionary tracks and 
spectral libraries. Following IWu et all (|2005[ l. we adopt 
"mod i = 0-05 in this paper. 

The BC03 SSP models include six initial metallicities, 
[Fe/H] = -2.2490, -1.6464, -0.6392, -0.3300, +0.0932 
(solar metallicity), and +0.5595. Spectra for other metal- 
licities can be obtained by linear interpolation of the ap- 
propriate spectra for any of these metallicities. For S312, 



whose metallicity and reddening values were published 
by other authors, the cluster age is the sole parameter to 
be estimated (for a given IMF and extinction law, which 
we assume to be universal). In Figure 3, we show the 
intrinsic SEDs of S312, the integrated SEDs of the best- 
fitting model, and the spectra of the best-fitting model 
superimposed the 3-sigma spectra. The best reduced % 2 
of 18.23 is achieved with an age of 9.5lJ;gg Gyr, which is 
in very good agreement w i th th e previous determination 
(10 Gyr) of lBrown et all (l200l by main-sequence pho- 
tometry. From Figure 3, we can see that, the 3-sigma 
spectra deviate clearly from best-fitting ones, i.e. these 
three spectra clearly separate. This indicates that we can 
constrain the age of S312 with errors of about 3 Gyr. 

3.5. Mass of S312 

Star cluster masses can be estimated by comparing 
the measured luminosity in V band with the theoreti- 
cal mass-to-light ratios. These ratios are a function of 
the cluster age and metallicity. BC03 calculated these 
ratios for a range of metallicities. Mass-to-light ratios 
for other metallicities can be obtained by linear interpo- 
lation of the appropriate ratios for any of these metallic- 
ities. For [Fe/H] = —0.53 and 9.5 Gyr, we derived the 
mass-to- light ratio in broad-band V to be 3.69 M Q /L . 
Based on its present luminosity, V — 16.13 ± 0.01 
mag and the extinction, E(B — V) = 0.10, its intrin- 
sic luminosity is Vb = 15.82 ± 0.01 mag [assuming the 
ICardelli et all (|1989f l Galactic reddening law; A v = 0.31 
mag]. From this, we determine the mass of S312 to 
be 9.8 ± 1.85 x 10 5 M Q . Co mparing with 037-B 327 
[A4 37-B327 ~ 8.5 x 10 6 M (IBarmbv et al.l l2002fl or 
A4 37-B327 ~ 3.0 ± 0.5 x 10 7 M ? dMa et al.ll2006f)l and 
Gl [M G i - (7 - 17) x 10 6 Mq (jMevlan et al.ll2001hl in 
M31 and w Cen [M^ C en ~ (2.9-5.1) x 1O 6 M ijMevlanl 
2002)] in the Milky Way, the most massive clusters in 
the Local Group, S312 is only a medium-mass globular 
cluster. 

3.6. Comparison with SSP models based on the 2000 
Padova evolutionary tracks 

As discussed in Section 3.1, BC03 provide two sets of 
SSP models based on the 1994 and 2000 Padova evolu- 
tionary tracks, respectively. In Section 3.4, we adopted 
the SSP models based on the 1994 Padova evolutionary 
tracks to constrain the age of S312. However, it is in- 
teresting to show what the best fit is based on the 2000 
Padova models. Using the 2000 Padova models, we fit- 
ted the intrinsic SEDs of S312 again, it is surprised that 
the fitting is very good with the best reduced \ 2 of 9.71, 
but, the age of 11.75+0 62 Gyr is older than 10 Gyr of 
iBrown et all (|2004f) by main-sequence photometry. In 
Figure 4, we plot the intrinsic SEDs of S312, the inte- 
grated SEDs of the best-fitting model and the spectra of 
the best-fitting model superimposed the 3-sigma spectra 
based on the 2000 Padova evolutionary tracks. This fig- 
ure also showed that, the 3-sigma spectra deviate from 
best-fitting ones much. As things stand now, without 
having other clusters to compare, there is a problem that 
an SED-fitting based on the older isochrones (the 1994 
Padova) constrains the age of S3 12 in agreement with 
the main-sequence turnoff fit, while the newer isochrones 
(the 2000 Padova) give an older age constraint of S312. 
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Fig. 2. — Comparison of the age scale between the Padova 1994 tracks and the 2006 vandenBerg isochrones used in Brown ct al. (2004). 
Dashed and solid lines are the Padova 1994 tracks and 2006 vandenBerg isochrone, respectively, for age 10 Gyr and the solar metallicity. 



These are therefore issues that one needs to keep in mind 
in the context of more recent input physical parameters 
of the 2000 Padova evo lutionary tra cks as indicated by 
IBruzual fc Charlotl (j2003h . In fact, iBruzual fe Charlotl 
(2003) have illustrated the influence of the stellar evo- 
lution prescription on the predicted photometric evolu- 
tion of an SSP for solar metallicity as an example in 
their Figure 2, which showed models computed using 
the 1994 Padova, the Geneva and th e 2000 Padova evo- 
lution ary tracks. From Figure 2 of IBruzual fc Charlotl 
(2003), we can see that, at late ages, the V — K color 
is significantly bluer in the 2000 Padova model than in 
the 1994 Padova model. The reason for this is that the 
giant-branch temperature in the 2000 Padova tracks has 
not been tested against observational calibrations (see 
IBruzual fc Charlotl 120031 and references therein). So, it 
is certain that an SED-fitting based on the 2000 Padova 
model gives an older age for S312 than the age we ob- 



tained from the 1994 IBruzual fc Charlotl (|2003h models. 

3.7. Overview of age constraint by photometry in 
different bands 

In section 3.4, we constrain the age of S312 based 
on the integrated photometric measurements in GALEX 
NUV, broad-band UBVR, 9 intermediate-band filters 
of BATC and infrared JHK S of 2MASS. These pho- 
tometries constitute the accurate SEDs of S3 12 from 
2267 - 20000A. In order to constrain the age of S312 
accurately, we use as many photometric data points cov- 
ering as large a wavelength range as possible. In this sec- 
tion, we investigate which passbands actually constrain 
the age. The results show that, except for NUV, a lack 
of any band does not result in great problems in con- 
straining the age of S312, the dispersions in the recov- 
ered age are smaller than 0.5 Gyr. The reason may be 
that, for S312, whose metallicity and reddening values 
were published by other authors, the cluster age is the 
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Fig. 3. — Plots of the best-fitting integrated SEDs of SSP models plotted on top of the intrinsic SEDs for S312. The photometry itself is 
shown by the symbols with error bars (vertical ones for uncertainties and horizontal ones for the approximate wavelength coverage of each 
filter). Open circles represent the calculated magnitude of the model SED for each filter. Yellow, red and blue indicate the spectra of the 
best-fitting model, plus and minus the 3<r spectra, respectively. 
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Fig. 4. — Graphs of the best-fit integrated SEDs of SSP models based on 2000 Padova models plotted on top of the intrinsic SEDs for 
S312. The symbols are the same as those in Figure 3. 
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sole parameter to be constrained. If we constrain the 
age and metallicity of S3 12 simulta neously, the i mpor - 
tance of the U or B will appear. As I Anders et alj (2004) 
pointed out that, the age deviations from the input val- 
ues for the combinations without the U or B bands are 
caused by an insufficiently accurate determination of the 
cluster metallicity, i.e. the U or B plays a major role 
in determining the metallicity. Missing U or £?-band in- 
formation leads to underestimates of metallicity, thereby 
causing ages to be adjusted improperly. The results also 
show that, a lack of NUV results in problems in con- 
straining the age of S312, the dispersion in the recovered 
age is 2.5 Gyr, which reflects that the NUV flux point 
shows promise of being useful determination of age. A 
lack of all bands bluer than V or of all bands redder 
than I results in great problems in constraining the age 
of S312, the dispersions in the recovered age are 4.25 or 
7.25 Gyr, respectively. 

4. SUMMARY AND DISCUSSION 

In this paper, we constrain the age of the M31 glob- 
ular cluster S312 by comparing its multicolor photom- 
etry with theoretical stellar population synthesis mod- 
els. This is both a check on previous results, and a 
check on our methodology of applying age constraints 
via our SEDs of stellar systems. Multicolor photomet- 
ric data are from GALEX NUV, broadband UBVR, 9 
intermediate-band filters and 2MASS JHK S , which con- 
stitute the SEDs covering 2267 - 20000A. Our result is in 
very goo d agreement with th e previous determination (10 
Gyr) of iBrown et alJ ((2004) using main-sequence pho- 
tometry We also estimate the mass of S312, which shows 
that it is a medium-mass globular cluster in M31. Since 
no other clusters can be available to confirm our SEDs- 
fitting method of constraining ages of globular clusters, 
a large sample of globular clusters, whose ages are de- 
rived based on our SEDs-fitting method, may boost con- 
fidence in the dataset and in our SEDs-fitting method. 



I Fan et aTl (|2006f ) obtain age constraints for 91 M31 glob- 
ular clusters, based on the same sources of photometric 
data (not including GALEX NUV) , the same theoretical 
stellar synthesis models and the same fitting methods. 
Figure 5 shows a plot of age error s (la spread) as a func- 
tion of age based on the data of iFan et alJ ((2006). It is 
clear that the absolute age errors increase significantly 
with increasing age, and derived ages with errors of < 3 
G yr are ava i lable for ages younger than 9 Gyr. 

IFan et all (|2006h argue for peaks in the age distribution 
at ~ 3 and 8 Gyr, in addition to the expected comple- 
ment of Milky Way-like "old" globular clusters. While 
their results have considerable scatter, they show a no- 
ticeable lack of young, metal-poor globular clusters, or 
old and very metal-rich globular clusters, which might 
be consistent with an age metallicity relationship. 

At the same time, the results of this paper also indi- 
cate that, despite having such good multicolor photom- 
etry, we still cannot get solid age constraints from SEDs 
fitting even with NUV to K s photometry. This has impli- 
cations for attempting to learn the ages of distant cluster 
systems. 



We are indebted to the referee for his/her thoughtful 
comments and insightful suggestions that improved this 
paper greatly. We thank S.-C. Rey for providing us with 
his data on S312 in GALEX NUV in advance of publi- 
cation. S.-C. Rey also kindly provides us the GALEX 
filter response curves. This work has been supported by 
the Chinese National Natural Science Foundation No. 
10473012, 10573020, 10633020, 10673012, and 10603006. 
This publication makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Pro- 
cessing and Analysis Center/California Institute of Tech- 
nology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 
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Table 1 

BATC Photometry of the M31 Globular Cluster S312. 

Filter A(A) FWHM(A) N a Magnitude 



e 


4925 


390 


11 


17.81(0.021) 


/ 


5270 


340 


12 


17.23(0.016) 


9 


5795 


310 


7 


16.40(0.012) 


li 


6075 


310 


5 


16.14(0.009) 


i 


6656 


480 


3 


15.54(0.007) 


3 


7057 


300 


12 


15.25(0.006) 


k 


7546 


330 


6 


14.89(0.006) 


o 


9182 


260 


18 


13.95(0.004) 


V 


9739 


270 


12 


13.78(0.005) 



a N is the number of images taken by the BATC telescope. 



Table 2 

nuv, broad-band and 2mass photometry for s312. 



Filter 


Magnitude Reference 


NUV 


21.881(0.123) Rev et al. ("20061 


U 


17.460(0.081 Battistini et al. (1987) 


B 


17.032(0.018) Reed. Harris & Harris (1994) 


V 


16.130(0.011) 


R 


15.560(0.014) 


J 


14.195(0.032) 2MASS 


H 


13.495(0.037) 


K a 


13.390(0.043) 



